Multiply Unstable Black Hole Critical Solutions 
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The gravitational collapse of a complex scalar field in the harmonic map is modeled in spherical 
symmetry. Previous work has shown that a change of stability of the attracting critical solution 
occurs in parameter space from the discretely self-similarity critical (DSS) solution originally found 
by Choptuik to the continuously self-similar (CSS) solution found by Hirschmann and Eardley. In 
the region of parameter space in which the DSS is the attractor, a family of initial data is found 
which finds the CSS as its critical solution despite the fact that it has more than one unstable mode. 
An explanation of this is proposed in analogy to families that find the DSS in the region where the 
CSS is the attractor. 
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I. INTRODUCTION 



Much progress has been made in understanding the 
nonhnear phenomena associated with the threshold of 
black hole formation first found by Choptuik He 
studied the collapse of a minimally coupled, real scalar 
field whose initial configuration is parameterized by some 
parameter p. For initial data with p less than some criti- 
cal value p* , the field energy explodes through the center 
of spherical symmetry and disperses. For p > p* , the 
configuration collapses to form a black hole. In the limit 
p ^ p*, a solution perfectly poised between collapse and 
dispersal, the critical solution, is reached. 

Naively one might expect to find a tremendous variety 
of critical solutions dependent on the form of the initial 
data used (eg. Gaussian pulses parameterized by an am- 
plitude p, sinusoidal pulses parameterized by a frequency 
p, etc). However, Choptuik found that the critical solu- 
tion for any interpolating family was precisely the same. 
The critical solution he found was said to be universal 
because of its apparent independence on the initial data. 

In addition to its universality, the critical solution 
found by Choptuik exhibits discrete self-similarity (DSS) 
such that the fields are periodic in log\t\ and logr via 
/(r, t) — /(exp(log r -I- A), exp(log \t\ + A)) for a univer- 
sal constant A « 3.44. Later, studies of axisymmetric 
vacuum gravity revealed similar critical behavior with 
a different echoing constant A « 0.6, showing A to be 
model-specific ||]. However, the nature of how the field 
equations select a specific value A was, and still is, mys- 
terious. 

Discovery of other critical solutions has since fol- 
lowed including various continuously self similar solu- 
tions (CSS) found in perfect fiuid collapse and other 
scalar collapse models whereby the fields obey f{r,t) = 

t^"/o(-rA) im. 

Here the issue of universality is examined. Univer- 
sality comes about because of the presence of only one 
unstable (relevant) mode about the critical solution. The 
single unstable mode sends nearby trajectories in phase 
space away from the critical solution. The trajectories 



either disperse or form black holes. This mode is then 
appropriately called the black hole mode. Though the 
critical solution is unstable, the process of tuning pro- 
gressively limits the infiuence of the mode, delaying its 
growth. Modulo this unstable mode, the critical solution 
is then an attractor and therefore called an intermediate 
attractor. 

When a critical solution has more than one unstable 
mode it ceases to be an intermediate attractor. Tuning a 
one-parameter family of initial data still tunes the black 
hole mode, but the trajectory will still generically be sent 
away from such a critical solution by the presence of the 
other unstable modes. 

The model studied here is the harmonic map (also 
called the nonlinear sigma model) with a free parameter 
At. The model represents a mapping of a self-gravitating 
complex scalar field onto a target space with constant 
curvature. This curvature, parameterized by —k, char- 
acterizes the internal space of the field. For the case 
K = the model is simply a free complex scalar field 
minimally coupled to gravity whose action we recover 
by setting k = in the action for the non-linear sigma 
model, Eq. (Al). In this case, the target space is the 
complex plane with zero curvature. For positive values 
of K, the non-linear sigma model maps into a hyperboloid. 
This model is equivalent to a real scalar field coupled to 
Brans-Dicke gravity studied in For negative k, the 
target space is the sphere, S*^. 

Because for k = the model is identical to the free 
complex scalar field, the attracting critical solution is 
known to be the DSS. Hirschmann and Eardley have 
shown the CSS to have a pair of conjugate unstable 
modes in addition to the black hole mode for n near 
zero However, their analysis shows that as k is in- 
creased above k « 0.075 a bifurcation occurs and the 
extra unstable modes become stable. Hence, in this re- 
gion, it should be an attracting critical solution, while 
below this range of k, the DSS is the critical solution. 
The evolutions of |^ confirm this change in stability and 
give evidence that the DSS is not the attractor above 
K ~ 0.1. The harmonic map then has both the CSS and 
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DSS as attracting critical solutions in distinct regions of 
parameter space. 

A remarkable family of initial data (called spiral initial 
data here) is presented which finds the CSS critical solu- 
tion in the region of parameter space for which the DSS 
is the demonstrated attractor. In this article, the reasons 
why this family can find a multiply unstable critical solu- 
tion are studied. It is argued that the spiral initial data 
is quite special in that its saturation of a charge bound 
disallows the growth of the extra unstable modes. 

This family is presented in Section || along with the 
critical solutions obtained. The family is then perturbed 
and the critical solution is seen to deviate from pure CSS. 
These results are discussed in Section [II , beginning with 
a discussion of the DSS in a region where it ceases to be 
an intermediate attractor. An explanation of the reasons 
that the CSS can be found when not an intermediate at- 
tractor follows. After the Conclusion, an Appendix pro- 
vides the equations of motion for the model. I also note 
that the results described below were generated using a 
modified version of Choptuik's adaptive mesh-refinement 
code for massless scalar collapse p|. 



F(r,0) = /(r)e^' 



(3) 



and where to is an arbitrary real constant. This data 
represents a generalized spiral in the complex plane and 
is called spiral initial data here. So that the fields are 
regular and compact, /(r) is constructed such that /(r — > 
0) = and /(rmax) = 0. The constant Tmax represents 
the size of the numerical grid. Also, for reasons that 
should become clear later, /(r) is constructed so that it 
varies much more slowly in r than the exponential term. 
Generally, / is either Gaussian or takes a step-function- 
like form 
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1 — tanh (r — r^^ 



(4) 



for arbitrary constants riow < '"high ■ 

Spiral initial data is remarkable because critical 
searches conducted with this data find the CSS critical 
solution for k « 0. Figure || displays the critical solu- 
tion obtained with various initial data for values of k 
for which the DSS is the attractor. These results show 
that the spiral data is quite special in the space of initial 
configurations. 



II. NUMERICAL RESULTS 



Because in spherical symmetry the gravitational field 
has no degree of freedom, initial data is completely de- 
termined by the specification of the scalar field and its 
time derivative at the initial time. Because the field 
F{r, t) is complex, we can decompose the field into its 
real and imaginary components ^(r,t) and (j){r,t)^ re- 
spectively, and specify their initial profiles '0(r, 0) and 
<j){r, 0). Setting the generalized time derivatives (see Eq. 
AlO) to zero initially via n^(r, 0) = n^(r, 0) = yields 
time symmetric data. Setting n^(r, 0) = il){r^Qy and 
n0(r, 0) = (/)(r, 0)' yields initial data which approximates 
an ingoing wave. For what follows, the choice of either 
of these does not affect the critical solution found. 

Generally the type of the initial data does not affect 
the obtained critical solution because of universality, and 
hence a common choice has been a Gaussian pulse in each 
of the components 



V'(r,0) 
0(r, 0) 



(1) 



where A^^ A^^R^^R^^d^^d^ are arbitrary real con- 
stants. However, instead of a decomposition into real 
and imaginary parts, the complex field can be expressed 
by a magnitude and phase 



F(r,0) = /(r,0) e"'('^^°) 



(2) 



for arbitrary real functions f{r,t) and h{r,t). 

The family which is of interest here is most easily ex- 
pressed in this form where the phase is linear in r 
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FIG. 1. Critical solutions for both Gaussian (dashed, Eq. 
|l]) and spiral (solid, Eq. ^ initial data for k = 0.1,0.0,0.5. 
The solution shown with open circles is that found by 
Hirschmann and Eardley with the assumption of continuous 
self-similarity [5]. 

Perturbations of this data are made according to 

V'('',0) = /(r) (A + AA) cos{ujr + Aip) 

0(^,0) =/(r) A sin(wr), (5) 

and the critical solutions obtained are shown in Figs. ^ 
and H The fi gures show the last subcritical solution ob- 
tained by a critical search at a time just before it de- 
cides to disperse. The graphs then represent the outgo- 
ing record of the collapse of the self-similar regime to- 
wards hir oo. Large r then represents early time, and 
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the graphs show that as the perturbation is increased, 
the self-similar pulse gradually develops a discrete oscil- 
lation. This development represents the tunneling of the 
solution away from the CSS and toward the DSS. 

These results make clear that changing the relative 
phase of the two fields or their relative amplitudes causes 
the critical solution to be attracted to the DSS. Pertur- 
bations of the relative frequency produces similar results. 



F(r,0) = /(r) exp (iwrP) 



(6) 



for p ^ 1 and still be considered 7r/2 radians out of phase. 
However, as shown in Fig. ^, for p 1 the CSS is not the 
critical solution. The initial data for these configurations 
are shown in Fig. ^. 
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FIG. 2. Perturbations of phase for the spiral initial data 
for K, — 0. Even for small changes in the relative phase of ip 
and (h, the critical solution is driven toward the DSS. 
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FIG. 4. Critical solutions (k — 0) obtained for the spiral 
data for various values of p [Eq. (H)]. The initial data for 
these critical solutions are shown in Figure ^ 
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FIG. 3. Perturbations of relative amplitude for the spiral 
initial data for k = 0. Again, the solution is driven toward 
the DSS. 

While Figures |^ and ^ show that only initial data com- 
pletely out of phase finds the CSS, the spiral data can be 
further perturbed via 



0.04 : p= 1.45 
0.02 V- 

— — 

-0.02 ^ 

0.04 E p=1.0 

-0.02 - 

0.04 ^ p-0.5 



dm/ dr 
V - - d(i/dr 

T 

--- 




-0.02 ^, 



10 



15 



20 



FIG. 5. Initial data for the critical solutions displayed in 
Figure ^ 

Before discussing these results, it is interesting to con- 
sider initial data which consists of the superposition of 
two different frequencies loi and uj2 

F{r,{)) ^ h{r)e'^'- + h{r)e'^^\ (7) 
Both the initial data and the critical solution are shown 
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in Fig. H This example of the superposition of two 
frequencies finds the CSS as shown in the figure, however, 
not all examples of this family do. Families with /i = /2 
and comparable frequencies resulted in critical solutions 
which appeared CSS. 
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FIG. 6. Critical solution obtained for initial data consisting 
of two frequencies ui — 2.1 and uj2 — 3.6 [Eq. (^)]. The top 
frame shows the initial data. The bottom frame displays the 
critical solution. 

As mentioned previously, initial data which is other- 
wise spiral but has /(r) varying on scales comparable to uj 
is driven away from the CSS critical solution in much the 
same way that the perturbations shown in Figures ^ and 
^ drive the solution away from the CSS. Also, many ex- 
amples of initial data of the form (|^) are similarly driven 
away from the CSS. Experimentally, the strongest indi- 
cator of initial data which will find the CSS is when the 
energy density is everywhere proportional to the charge 
density. That this proportionality indicates the spec ial- 
ness of spiral initial data is discussed in Section 111 B| . 



III. DISCUSSION 

Surprisingly enough, in this case a discussion of the 
DSS critical solution in the region k > 0.1 is simpler than 
the discussion of the CSS. In this region of parameter 
space, it is the DSS which has multiple unstable modes, 
and it is relatively easy to understand the families of 
initial data which find the DSS. Hence, the discussion 
of these special families is presented first, followed by a 
discussion of the specialness of the spiral data. 



A. The DSS 

In the harmonic map model, for k, ^ O.l the CSS is 
the attracting critical solution. Numerical evolutions of 
various families of initial data generically find the CSS, 
and do not find the DSS. Gundlach showed that the DSS 
has only one unstable mode for k — 0, but the stability 
analysis of the DSS has not been done for general k . 
However, evolutions of this model, as well as results in 
the equivalent region of the model in j^] , clearly indicate 
that for K > 0.1 the DSS has more than one unstable 
mode. 

However, there are non-generic families that do find the 
DSS in this regime despite the presence of these extra 
unstable modes. Consideration of this phenomenon is 
helpful in understanding the CSS occurring for k = 0. 

One description of initial data that finds the DSS in 
this regime is mentioned in ||^ . In that work it was found 
that when one component of the field was initially vanish- 
ing, it remained zero. The model here, being equivalent 
to that one for k > 0, retains this feature as shown in 
the equations of motion for the two components of the 
scalar field in Equation (A8). Because the CSS is neces- 
sarily complex (it has charge) , initial data with one field 
initially vanishing is unable to find the CSS as its critical 
solution. Thus, families of initial data of the form 



F{r) = f{r) 
iP{r) = f{r) 
<f>{r) = 



(8) 



for arbitrary /(r) will only find the DSS. 

However, a more general set of families can be found 
with this principle in mind. Consider initial data of the 
form 



JC 



f{r) 



Fir) 

V'(r)=cos(C)/(r) 
0(r) = sin (C)/(r) 



(9) 



for arbitrary constant C. This data corresponds to a 
global rotation of Eq. (||) by an angle C in the complex 
plane. Because the Lagrangian is invariant with respect 
to this rotation, the critical solution must be the same as 
the initial data described by Equation (H). 

A more physical understanding of this can be gained 
by examining the issue of charge. For both sets of initial 
data (HH), the charge is zero. In fact, all components of 
the current density. Equation (A3), vanish 



= 0. 



(10) 



The divergence of this current is zero, so the current den- 
sity will not grow if initially vanishing. In other words, 
the system with no charge is in a symmetric state with 
respect to charge, and this symmetry would have to be 
broken were the charge to become positive or negative. 
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Because the CSS can have either positive or negative 
charge, these sets of initial data are precisely balanced 
between trajectories that would take them to the CSS 
with positive charge and those that would take them to 
the negatively charged CSS (since the model is indepen- 
dent of which field is considered the imaginary compo- 
nent and which the real component of the complex scalar 
field, the sign of the charge of the CSS is arbitrary). It is 
this balance that enables them to see the DSS as a critical 
solution with only one unstable mode when in fact it has 
more than one. The initial data has already tuned one of 
the extra unstable modes (or an unstable conjugate pair 
of modes). 

With the knowledge that the extra unstable mode cor- 
responds to a charged mode, a two-parameter search 
can now be conducted. Determination of an appropri- 
ately parameterized initial data family is somewhat more 
subtle than that for the one parameter data. With a 
one parameter search, a parameterization needs to be 
smooth and monotonic in the mass of the initial data 
near the critical point because the excitation of the black 
hole mode is characterized by the mass contained. Here, 
that mode must also be tuned, but the unstable charged 
modes must be tuned as well. Hence, the second param- 
eter must be locally monotonic in charge near the critical 
point. 

A two parameter search with the data 



^(r) =pf{R^) 

4>{r) f{R^-6 + 10p) 



(11) 



finds the DSS. Here, f{R) is a Gaussian pulse centered 
on some radius R, and R^ and D are arbitrary constants. 
The parameters p and S are used to tune both the mass 
and the charge of the initial data. For 6 fixed, p effectively 
tunes the initial energy content. With p fixed and with 
6 ~ lOp, increasing d increases the charge from some 
negative value up to some positive value (the constant 10 
is chosen arbitrarily). Setting S = lOp turns the family 
back into one with zero charge. 

With this data, some value for S near the value lOp is 
chosen, and the black hole critical solution is bracketed 
sufficiently closely so that the sign of the initial charge 
of the critical solution is known. A different 6 is chosen 
so that the sign of the charge of the critical solution is 
found to be the opposite sign. These two values of 5 then 
bracket a critical value 6* for which the critical solution 
{p — p*) has zero charge. In this manner, the DSS is 
found. 

The difference between this two-parameter tuning and 
using initial data with zero charge is simply the order in 
which the tuning occurs. In the latter, the initial data is 
already tuned to have zero charge. The remaining task is 
then to tune the mass. However in the former, the charge 
is being tuned first to arrive at a one-parameter family 
that in general has charge. It is only at the critical value 



of p that this one-parameter family (found from the 5 
search) that the solution has zero charge. 

In Figure ^ a series of critical solutions for n = 1 with 
initial data of the form (|ll|) are shown. For each case, 
\p^p* I is of order machine precision. The critical solution 
gets closer to the DSS in the limit [5-5*) /5* 0. In- 
terestingly, while the critical solutions of the perturbed 
spiral data in Figure ^ appeared initially to be headed 
toward the CSS only to be funneled to the DSS, these 
solutions do the opposite. They begin initially as DSS 
solutions but eventually head to the CSS. As they get 
closer to the critical solution (5*, the DSS lasts for a pro- 
gressively longer time. This similarity is consistent with 
the spiral being tuned to the unstable modes of the CSS. 
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FIG. 7. Critical solutions with n = 1 for various levels of 
charge tuning (5 — 5*) /5*. The initial data is of the form 
found in Eq (|ll|). The field 2m /r is shown for {5-5") /S* 
equal to 0.1,0.01,0.001,0.0001. As 5 approaches 5* the solu- 
tion appears discretely self-similar longer but eventually be- 
comes the CSS. 



B. The CSS 

Initial data of the form Equation (^) will always find 
the DSS, because, having no charge, a direction must 
be picked toward either positive or negative charge to 
break the symmetry. In the case of finding the multiply- 
unstable CSS, the spiral data maximizes the charge of 
the initial data for a given energy, and so it appears that 
a symmetry must be broken to disperse all the charge 
and arrive at the DSS. 

To find a bound on the charge, consider the norm of 



the vector 



papapc^ ± 



at the initial 



time for time-symmetric initial data. Because at the ini- 
tial time all the time components vanish, the norm is 
positive definite. This is similar to a trick employed by 
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Belavin and Polyakov g and described in Rajaraman | 
Computing 



(l — k\F\^) 



papa pc,fi _|_ 



)^PYp^^p*,t^ > 



> (12) 



(13) 



This inequality ( [l3| ) can be expressed in terms of the 
energy density for time-symmetric initial data 



Tqo — 



a^F„F*^^' 



87r(l-K|F|2)^ 
The bound is then 

2,2tt\F? 



Too > a'3^,f 



(1-^|FP) 
which, for k = 0, is simply 



(14) 



(15) 



(16) 



Hence there is an upper limit to the square of the charge 
density for a given energy density. 

Letting the initial data be of the general form (H), the 
condition to saturate this bound for k = is 



which implies 



(/') = 0. 



(17) 



(18) 



Saturation therefore occurs when all energy occurs in the 
phase rotation, h' . Physically, this is apparent by looking 
at the behavior of F in the complex plane. Because /' 
vanishes, as F is traced out for various r, the magnitude / 
does not change; only the phase is changing so the path 
is a circle on the plane. This tracing then maximizes 
the area covered for a given energy. The area covered is 
proportional to the charge, so the charge is maximized. 

This analysis applies only for At = with time- 
symmetric initial data, though presumably similar argu- 
ments would hold for the generalization to other values 
of K and ingoing initial data. 

At this point, it is interesting to compare the family (|^) 
to the spiral initial data. These two families are, in some 
sense, complementary. Pick some r, and imagine that 
point in the complex plane. To construct data of the 
form (^, determine the field values for all other r by 
requiring these points to fall on the radial line between 
this initial point and the origin. In this fashion, the initial 
data will have some global phase, constant in r, equal to 
some value C. However beginning once again from that 



initial point in the plane, the restriction of Equation ( [T^ ) 
says that to construct spiral initial data, as r is increased 
the curve in the plane must lie everywhere perpendicular 
to the radial direction. 

Another condition restricting the spiral data appears 
to be that the charge density must be independent of r. 
Just as seen with families that find the multiply unsta- 
ble DSS where the charge density must be everywhere 
zero independent of r, here the charge density must be 
everywhere a maximum and independent of r. 

In order for the charge density to be independent of r, 
the rate at which the data covers the plane, lu — h' , must 
be constant in r. Initial data which would otherwise be 
spiral but with h' ^ are shown in Figure ^ with their 
respective critical solutions shown in Figure ^. These 
results indicate that only for p = 1 in Equation (^ is the 
CSS found. 

Another perspective on the constraints of the spiral 
data is afforded by examining the ratio of the charge 
density to the energy density 



2a^f^h' 



F^,F*'t^ 



ifr + p {h'Y 



2a2 
h' ■ 



(19) 



That this ratio is approximately independent of r for the 
conditions /' = and h" = appears to indicate that 
tuning the mass of the initial data also tunes the charge. 
In other words, with the charge to mass ratio everywhere 
the same, the critical search cannot find a solution which 
disperses all the charge. 

Construction of smooth, compact, and regular initial 
data consistent with the restrictions 



/'(r,0) = /i"(r,0) = 



(20) 



is quite difficult. Regularity at the origin requires either 
/(O,i) = 0or f{0,t) = = h'iO,t). 

Satisfying the former along with strict observance 
of ( [20| ) leads to the trivial solution F(r,0) = 0. Instead, 
as mentioned earlier, initial data is used where / vanishes 
at the origin, but is then "turned on" at some larger r. 
The condition /' = is then not satisfied everywhere, 
but for /' < w the CSS is still found. 

Satisfying the latter condition along with h" = at the 
initial time also leads to a trivial solution F{r, 0) = C, for 
some complex constant C. Again, to can be "turned on" 
at some larger r, but this can only be done in a smooth 
way if h"{r, 0) is not everywhere zero. 

The difficulty in constructing non-trivial, regular, ini- 
tial data which is strictly spiral has hampered the analy- 
sis. While initial data of the form (^ has the symmetry 
(j){r) — tan (C) ip{r) which holds at all times, it is not 
clear if there is such a preserved symmetry here. The 
approximate symmetry 0(r) = tan (wr) ipir) holds at the 
initial time for the spiral data but does not appear to 
hold throughout the evolution. Also, the bound (|8|) is 
shown only for time-symmetric initial data but ingoing 
spiral data also finds the CSS. However, the evolutions 
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consistently show that initial data which has the ratio 
of charge density to energy density independent of r will 
find the CSS as its critical solution. 



IV. CONCLUSION 

A special family of initial data is described which finds 
the CSS critical solution in a region of parameter space 
where generic initial data finds the DSS. The specialness 
of this family is discussed. A bound on the charge den- 
sity for time-symmetric initial data is found, which this 
family saturates. Because the spiral data maximizes the 
charge density and because this charge density is inde- 
pendent of r, it is argued that the extra unstable modes 
cannot grow because their growth would pick a direction 
in which to decrease the charge. The inability of the 
extra modes to grow then indicates that the spiral data 
"see" only the black hole mode around the CSS. 

Because both the CSS and DSS can both be found 
where they have multiple unstable modes, it seems pos- 
sible that other suitably tuned families in other models 
might find other, previously unknown multiply unstable 
critical solutions. 
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APPENDIX: EQUATIONS OF MOTION 

The action S for the model under study here is 

2|Vi^|2 \ 



S = / (fx^/^ R - 



(Al) 



{l-K\F\^Y) 

defined in terms of a complex scalar field F(r, t) 

F{r,t) =ij{r,t)+i(l){r,t), (A2) 

its complex conjugate F*{r,t), and a dimcnsionlcss pa- 
rameter K. The operator V represents the covariant 
derivative. The action (Al) is invariant with respect to 
global rotations of F, and thus has a conserved current 



i{FF\^-F*F^^) 
{l-K\F\^f 



(A3) 



In component form where Latin indices run over 1 and 2 
for the real and imaginary components, this current is 



{l-K\F\^) 

The field equations are then 



-2kF* 

OF 



(A4) 

(A5) 
(A6) 



where G^^ is the usual Einstein tensor and the stress 
energy takes the form 



47r(l-K(?A2^02))2 ■ 



(A7) 



In terms of the real and imaginary parts of F, the wave 
equation becomes 



U(h = - 



1-k(V'2-H02) 



(A8) 



1 - K (^^2 -I- 02) 

We work in spherically symmetry with the metric 

ds^ ^ ~a^{r, t)dt^ + a2(r, t)dr^ + r^dil^. (A9) 

We introduce the following auxiliary variables in order to 
cast the field equations in first order in time form 



Hw, = -Tp 



H. 



a 
a 



(AlO) 



where overdots and primes denote derivatives with re- 
spect to t and r, respectively. The two second order wave 
equations become four, first order equations 



(All) 



2aK 



n - ^ + - n^n^) 



H, 



a(l - niip^ + cj)'^)) 



r a 



2aK 



(A12) 
(A13) 

($2 _ n2 _ $2 + n2 ) + 2^ ($^$^, _ n^H, 



a(l-K(l/'2 + (^2)) 



(AM) 



The fields ijj and 4> are maintained at each time step by 
spatially integrating their respective spatial derivatives 
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(j){r, t) 



$0(f,t) df. 



The Hamiltonian constraint is 



a 
a 



2r 



2(1-k(Vj2+02))2 



(A15) 



(A16) 



(A17) 



The nature of polar shcing and the radial gauge yields 
the constraint on the lapse function a 







a 
a 



a 
a 



1 



(A18) 



Finally, combination of an evolution equation and a mo- 
mentum constraint yields an evolution equation for a 



(1 - K (V'^ + (/.2))2 ■ 

Regularity of the solution at the origin demands 

a{0,ty = a(0,<)' = 
$^(0,i) = $0(O,i) = 0, 



(A19) 



(A20) 
(A21) 



and local flatness there is enforced by 
a(0,<) = 1. 



(A22) 



We have the freedom to pick a condition on a on each 
time slice which corresponds to a global change in the 
labeling of slices. We impose the condition on a at the 
large radius boundary of the grid rmax 



1 



(A23) 



As long as no radiation is escaping from the grid, this 
condition implies that coordinate time t corresponds to 
proper time for an observer at r = cx). 

The nature of being restricted to a finite grid imposes 
the need for an artificial boundary condition on the mat- 
ter fields there. Since our spacetime is asymptotically 
flat, we impose an approximate outgoing radiation con- 
dition on the matter fields. The flat space wave equation 
for some general scalar field 0(r, t) in spherical symmetry 



□e : 
has the solution 

re = 



(re),, = (re). 



f{r + t)+ g{r - t) 



(A24) 



for two general functions / (in-going component) and g 
(out-going component). To eliminate the in-going com- 
ponent at the outer boundary, we enforce at the bound- 
ary the condition 



9 (re) _ a (re) 



dr 



dt 



with the equations 



ne + n'e + ^ = o 

$e + + — = 0, 



(A26) 

(A27) 
(A28) 



which effectively limits reflection off the outer boundary. 

The metric ( |A9| ) corresponds to a "dynamical" 
Schwarzschild metric, allowing the association 



.„,o^(i-^)"", 



(A29) 



where the field m{f,t) represents the mass aspect func- 
tion measuring the amount of mass contained within a 
shell of radi us f c entered about the origin at coordinate 
time i. Eq. (|A29|) leads to 



/ , r / 1 



(A30) 



so that from the field a{r,t), we can obtain the mass 
aspect function. 

An advantage of polar slicing (and most slicing condi- 
tions used in numerical relativity) is that it avoids singu- 
larities. Further, because of our use of polar-areal coordi- 
nates, we cannot observe the formation of a true horizon 
which is a coordinate singularity in our coordinates. In- 
stead, by monitoring m(r, t)/r we can observe the forma- 
tion of a horizon and hence a black hole when m/r 1/2. 
Where m/r approaches 1/2, we can determine the mass 
of the black hole forming by simply halving the radius 
where the horizon forms. 
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